Astronomy & Astrophysics manuscript no. nonlinrev February 2, 2008 

(DOI: will be inserted by hand later) 



On non-linear hydro dynamic instability and enhanced transport 

in differentially rotating flows 

Denis Richard 1 ' 2,3 

1 NASA Ames Research Center, MS 245-3, Moffett Field, California 94035 

2 LUFH, Observatoire de Paris, F-92190 Meudon Cedex 

3 GIT-SPEC, Commissariat a l'Energie Atomique, l'Orme des merisiers, F-91190 Gif-sur-Yvette 
Received ; Accepted June, 18 2003 

Abstract. In this paper we argue that differential rotation can possibly sustain hydrodynamic turbulence in 
the absence of magnetic field. We explain why the non-linearities of the hydrodynamic equations (i.e. turbulent 
diffusion) should not be neglected, either as a simplifying approximation or based on boundary counditions (?). 
The consequences of lifting this hypothesis are studied for the flow stability and the enhanced turbulent transport. 
We develop a simple general model for the energetics of turbulent fluctuations in differentially rotating flows. By 
taking into account the non-linearities of the equations of motions, we give constraints on the mean flow properties 
for the possible development of shear instability. The results from recent laboratory experiments on rotating flows 
show - in agreement with the model - that the pertinent parameter for stability appears to be the Rossby 
number Ro. The laboratory experiments seem to be compatible with Ro < — 1/2 and Ro > 1 in the inviscid or 
high rotation rates limit. Our results, taken in the inviscid limit, are coherent with the classical linear stability 
analysis, in the sense that the critical perturbation equals zero on the marginal linear stability curve. We also 
propose a prescription for turbulent viscosity which generalize the /3-prescription derived in ? 
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1. Introduction 

Differentially rotating flows are present in a wide variety 
of Astrophysical systems, including stellar interiors, ac- 
cretion disks, or gaseous planets. Analytical studies teach 
us that differentially rotating flows are unstable, accord- 
ing to linear stability theory, whenever their angular mo- 
mentum is decreasing outward (?). These flows have been 
extensively studied but they do not exist in astrophysi- 
cal context. Nevertheless, linear stability does not mean 
stability. This is because turbulence is genuinely a non- 
linear phenomenon. Even when the primary instability is 
linear in nature, the self-sustenance of fully-developed tur- 
bulence requires for the non-linearities to enter the game 
at least - but not only - for the saturation of amplitudes. 
In that sense, linear stability analysis is limited to predict, 
for a restricted class of flows, their instability and only the 
very early development of the bifurcated solution. It also 
means that linear analysis predicts instabilities but does 
not predict unconditional stability. Numerical models fail 
to maintain hydrodynamical turbulence in keplerian flows 
(?). Based on past laboratory experiments from the 1930's 
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(??), differentially rotating flows can become unstable for 
Reynolds numbers of order 10,000, a figure that to our 
knowledge has not been reached in published numerical 
simulations of keplerian flows. 

? have shown, using laboratory experiments published 
in the 1930's, that differentially rotating flows exhibit 
finite-amplitude instabilities in the case where both an- 
gular velocity and angular momentum are increasing out- 
ward, a class of flows relevant for stellar interiors mod- 
els. Recent laboratory studies (?), showed that the class 
of flows experiencing shear instabilities is much wider, as 
demonstrated by experimental results. We aim in this pa- 
per to find some hints about the physical mechanisms in- 
volved in these instabilities. 



2. Background and motivations 

In this section we present the basic equations for the evo- 
lution of turbulent fluctuations. We discuss the relevant 
approximations that can be done and their implications 
regarding the properties of the turbulence and the mean 
flow. 
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2.1. Turbulent transport equations 

We consider the equation of motion for a viscous incom- 
pressible flow, 



dtU + (U.V)U = --VP + vV 2 U 

P 



(1) 



and we decompose the flow as its mean velocity and pres- 
sure fields and their fluctuations, namely 



U = u 



u 



P = P + P 



(2) 



where U and P are the total velocity and pressure fields, u 
and p are the time averaged velocity and pressure and u' 
and p are their time dependent fluctuations. Multiplying 
the equations of motion by the fluctuation field and aver- 
aging over time, considering that 



u = 



P =0, 



(3) 



(where x is the time average of x) we find the evolution 
equations for the turbulent fluctuations kinetic energy (i.e. 
the diagonal components of the Reynolds stress tensor) 
reported in appendix 1X1 Note that viscous diffusion has 
been neglected, but viscous dissipation has been kept. We 
consider circular motion = = 0), with axial and 
azimuthal symmetries (d z =8^ — 0) for the mean fields 
and define — = fi. The time averages of axial and az- 
imuthal derivatives of fluctuating terms is neglected (i.e. 
net turbulent transport is only radial). 



1„ — 



2fl u'j.u^ -d r (ru' r ) + 

+ - (j)'d r u' r - drp'u'^j 



2^ , U r u <p 



^(d r u' r r + ^(d r u'^ + ljd r u' z )^, (4) 



2 dtU * 



2 z 



d r ftr 



U rU(j) 



^(^(^;) 2 + KA) 2 + ^KA) 2 )' ( 5 ) 



-fd r {ru' r u' z ) - -{d r u' z ) 2 . 



(6) 



This set of equations can be simplified by making the 
assumption that the turbulence is close to isotropy. This 
allows us to neglect the time averaged pressure fluctua- 
tions terms (see for example ?). We further assume that 
the characteristic spatial variation scale of the fluctuating 
components is smaller than the local radius (i.e terms of 
order u^^u^jr can be neglected compared to the ones of 
order drU^^Uf,). The above equations then become, 



1.— 



2 * 

u- 

2 z 



20, u' r u^ — u' r d r (u' r ) 

- v \jd^+ ^W<A 

OMr 2 



u'ru'tt, - u'd r {u' r u') - -{d r u'f 



-u z d r {u' r u' z ) - -(d r u' z ) 2 . 



The total turbulent kinetic energy, defined as 
1 
2 



then obeys the following equation 

d t k = —u'rU.rdrVL + u i d r u r u i — v{d i u i ) 2 . 



(7) 

(8) 
(9) 

(10) 
(11) 



The first term of the rhs in equations Q) and © rep- 
resents the production of turbulence. We will also refer to 
them as "coupling" terms, as they couple the turbulent 
fluctuations to the mean flow. The second term of the rhs 
in the same equations and the first one of equation 
describe the diffusion due to turbulence. The last term of 
equations 0-© is the viscous dissipation. 

2.2. About linearized equations 

In this section we will be considering the effect of ne- 
glecting the non-linear terms (i.e. turbulent diffusion) in 
Equations Q-injl and There are several reasons why 
one would neglect the non-linearities of the hydrodynamic 
equations. The first one is to simplify the system by ne- 
glecting terms that are difficult to treat in an analytical 
stability analysis. Other reasons might include the belief 
that the physical system is such that the non-linearities 
actually have a negligible effect on its local or global evo- 
lution. Neglecting a priori such terms based on their local 
effect is rather difficult, considering that little is know on 
their actual amplitude which does depend closely on the 
intrinsic nature of the turbulence. Another option is to 
consider spatially averaged equations and make some hy- 
pothesis on the boundary conditions of the systems. For 
accretion disks, it has been argue that the average over 
the whole flow or part of the flow of the turbulent diffu- 
sion terms eventually vanishes. Making this hypothesis is 
equivalent to saying that there is no energy flux through 
the boundaries of the system. This is valid for an isolated 
system, but it seems unlikely for a disk, as it would not al- 
low for energy transfer with its surroundings, in particular 
the central star. It does not seem realistic just by consid- 
ering the mass flux between the two objects. It is likely 
that such energy transfer occurs through the boundary 
layer between the star and the disk. Imposing that condi- 
tion also formally implies that the turbulent stress tensor 
vanishes on all boundaries, meaning that there is no tur- 
bulent viscosity at the edges of the disk. Which is more, 
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the radial velocity also has to vanish at the boundaries. 
We could also consider the properties of the boundaries 
beyond which the system can not be considered as a con- 
tinuous medium, hence cease to be described accurately 
by the equations of hydrodynamics. The description of the 
interactions between the fluid component of the disk and 
its outer non-continuous part is not a straightforward task. 
Our last argument on this matter is to point that these 
turbulent diffusion terms are the ones that are modeled 
by the anomalous turbulent viscosity. Neglecting them is 
equivalent to making the assumption that there is no tur- 
bulent viscosity, no enhanced transport. It is then an in- 
evitable result, that starting with such a truncated set of 
equations, one reach the conclusion that there is no such 
transport. 

Considering only the linear terms of the turbulent 
transport equations will give us insights -by definition- 
only on the linear stability of the flow. The linearized set 
of equations is : 



Table 1. Correlation coefficients properties (Inviscid case) 
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3. Considerations on the non-linear stability 

In the following sections, will focus on the case of linearly 
stable flows (i.e. where d r Vlr 2 > 0) and consider the energy 
equations including the turbulent diffusion. 



3.1. Inviscid flow 

We first consider the case where dissipation due to molec- 
ular viscosity can be neglected. We define u the charac- 
teristic turbulent velocity (e.g. the root-mean-square ve- 
locity) and A the length scale characteristic of the spatial 
variations of the velocity fluctuations, and we pose 



— OtUj, 

2 * 



2f2 u'j.u'^ 



v (W) 2 + \(d r u^ + \{d r u z )^j , (12) 

(13) 



d t k = —u'rU^rdril — K<^ M i) 



'\2 



(14) 
(15) 



From Eq. I|15[) we conclude that a stationary turbulent 
state (dtk = 0) can exist only if the following condition 
-depending on the sign of the angular velocity gradient- 
is satisfied 



d r n > o 



u' r u, < , 



d r tt < => u'rU^ > . 



(16) 



For Rayleigh unstable flows (where d r flr 2 < and 
therefore dtt < 0), eqs. (O, (H2J) and JTH) show that the 
linear coupling terms with the mean flow are always a 
source for turbulent fluctuations, a property that reflects 
the global linear instability of this class of flows. The non- 
linearities are needed there only for the saturation of am- 
plitudes of the radial and azimuthal velocities along with 
the redistribution of energy towards axial motions of the 
new bifurcated flow. In the case of stable angular momen- 
tum stratification (<9 r £lr 2 > 0), the first order coupling 
terms have opposite signs in Eqs JJJ and (JBJ). This means 
that one of them is an energy sink for one of the com- 
ponent of the velocity fluctuations (?). It then looks im- 
possible to allow for the growth or to maintain turbulence 
as long as the angular momentum is increasing outward. 
This actually reflects the linear stability of such flows. 



d t u' 2 = 2n CL u 2 



a 



— d r VLr 2 , 
2 dtU 4> = 



r<p 



T' 

: + c: 



A ' 



1 — 



u 3 



d t k = -rdM C, 



T(f> 



fin U 

c y 



(17) 
(18) 
(19) 
(20) 



We have introduced the correlation coefficients C™ z , 
C$,Cl l , C l r4> and C n = C? + C$ + C?. By imposing simple 
constraints, we can deduce the signs of these parameters, 
which depend on the sign of the angular velocity gradient. 
In the case of a negative gradient for example, we know 
from Eq. l|15fl that the correlation product u' r u^ has to 
be positive, which translates into C l r j > 0. It also implies 
that a stationary state can be reach only if C™ is positive 
and is negative. Note that C™ must be zero in that 
case because we neglected the viscous dissipation. These 
properties for all introduced coefficients are summarized 
in Table . From Eq. (|20|l we derive the amplitude of the 
fluctuations vorticity : 



A 



a 1 



(21) 



This relation express that the vorticity extracted from 
the background flow by the fluctuations is proportional to 
the local shear. 

We have seen that there always exist a sink term for 
one of the components of the energy fluctuations. A nec- 
essary condition for the existence of self-sustained turbu- 
lence is then that the non-linear terms overcome this neg- 
ative production term (for the azimuthal component when 
d r fl < 0, for the radial one when d r il < 0). Namely, from 
Eqs JTIJ} and (HHJ 
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d r n > o 
d r n < o 

Which reduces to 



C?y >-2flC^ u 2 



Cij, u 2 



(22) 
(23) 



d r n > o 

dM < 



A C™ 

it d r £lr 2 



(24) 
(25) 



Finally, using Ea. (|21J) . we obtain the following rela- 
tions, 



d r n > 
d r n < 



fil fit 

Cr0 a „ . C% <9 r f2r 2 
— -ro r \l ^ — — 



(26) 
(27) 



Introducing the dimensionless Rossby number, Ro = 
rd r £l/2Q, and using the property C n = Cfi+CV; we obtain 
the necessary conditions for instability : 



d r Q > 
d r Q < 



Ro>-— = Ro c 



(Ro c > 0) , (28) 
(i?o c <0) . (29) 



Where we have introduced the parameter Ro c . A sim- 
ilar result can be derived by simply posing dt < u i >= 
and reducing equations fTTjl and lfH3|) . which leads to 
Ro = —C n /C™. The quantity Ro c can be seen as con- 
stant by assuming that the ratio C"/C™ is independent 
from the mean flow. This ratio quantify the redistribution 
of the energy extracted for the mean flow between the ve- 
locity components, by the non-linearities. Note that the 
coefficients C™, C^,C™ and C n themself most likely de- 
pend on the mean flow (see ?). Within this picture, Ro c 
is a constant critical stability parameter, and we expect 
the flow to become turbulent for either Ro > Ro c (when 
d r il > 0) or Ro < Ro c (when d r Vt > 0). ? already noted 
that the Rossby number should approach a constant value 
when a linearly stable rotating shear flow becomes turbu- 
lent, by arguing that " it is conceivable and even likely that 
an asymptotic state can be reached with a constant ratio of 
transfer from the transverse motion to total energy trans- 
fer from the mean flow" . The two actors in the balance 
are the gradient of angular velocity, shearing the velocity 
fluctuations, creating small scale vorticity, feeding the tur- 
bulent cascade, and the gradient of angular momentum, 
stabilizing the flow by damping the fluctuations through a 
" spring" mechanism. The energy extracted from the back- 
ground flow into small scale vorticity must reach a critical 



Table 2. Correlation coefficients properties (Viscous case) 
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value in order to overcome the constraint introduced by 
the stable gradient of angular momentum. The Rossby 
number, as a measure of the ratio between the shear and 
the angular momemtum gradient, appears a posteriori as 
a natural control parameter of the flow stability. 

3.2. Viscous flow 

Following the same path as in the inviscid case, we now 
add the viscous constraints to equations (|f 7|l - l|20() : 



f 



d t u r 2 = 2n c l rA u 2 + c; 1 ^ - vc?^, 



i — 



J r<$> 

d r nr 2 
r 

,3 



C% u 2 



2 t z ~ X A 2 ' 



3 2 

d t k = -rd r n C l r4> u 2 + C n y . - vC v ^ 



(30) 

(31) 
(32) 
(33) 



where C v = C v r + C% + C v z . The properties of the cor- 
relation factors are reported in Table (J2J). The extraction 
of vorticity from the mean flow now becomes : 



- - 9l± R O-L _ C " 

A ~ C" r A 2 C" ' 

and the relations l|24|l and l|25|l translate into 



(34) 



d r n > o 

d r Q < 



n\2 ' 



CJ} A 



A ' C™ r C"A 2 



(35) 
(36) 



and after some rearrangement we finally obtain, 

d r n > : Ro > Ro c (l - l/Re+) , 
d r n < : Ro < Ro c (l + l/i?e~) , 

where we have introduced the following quantities, 

2fl\ 2 



(37) 
(38) 



Rer = 



v 

2nx 2 



v 



{C u JC l r4> - C v C^IC n C l r ^)- x < 0, (39) 
(C2/CL - C v C n JC n C')- x > 0. (40) 
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Fig. 1. Stability diagram of Couette- Taylor flow in the (0 o ,0j) plane where o and O, are the angular velocities (in 
Hz) of the outer and inner cylinder respectively ; solid line : solid body rotation ; dotted lines : viscous stability criteria 
; dashed lines : inviscid stability criteria (see text for detail). On the right (resp. left) of the solid body rotation curve, 
the angular velocity gradient is positive (resp. negative). The flows within the dark surface are Rayleigh unstable. The 
dots are the experimental instability boundaries within which the flows exhibit self-sustained turbulence ; The lined 
surfaces stand for the predicted non-linearly instability. Experimental results are from ? 



The correction due to viscous dissipation lies in the 
1/Re + and l/Re~ terms. The necessary condition for in- 
stability reduces to the results found for the inviscid case, 
if v is set to zero. For high rotation rates, the critical values 
of the Rossby number will eventually tend to the inviscid 
values. 



4. Turbulent transport and viscosity 

Radial turbulent transport of angular momentum is quan- 
tified by the second-order correlation product u' r u'^. The 
time evolution for this quantity after some rearrangement 
and spatial averaging over the azimuthal direction <j>, and 
the introduction of correlation coefficients, is given by : 



Considering all of the above, we expect that there ex- 
ists a critical Rossby number for each class of flows - 
namely flows with <9 r O > and flows with <9 r O < 0. This 
result can be verified by using experimental results on dif- 
ferentially rotating flows. This can be done by examin- 
ing the value of the Rossby number at the onset of shear 
turbulence observed in laboratory Couette- Taylor experi- 
ment. The Couette- Taylor experiment consist of two coax- 
ial cylinder between which the fluid is sheared, resulting 
in a differentially rotating azimuthal flow. Experimental 
results on stability (?) can actually be fitted with the val- 
ues Ro c = 1 for d r Q, > 0, Ro c = —1/2, for <9 r O > 0, 
Re + — —0/3 and Re" = 20. This is consistent with the 
picture of a constant critical Rossby number with a vis- 
cous correction vanishing for high rotation rates or low 
viscosity. 



- -c^u 2 rd r n + c^- - v c%- 2 + n r0 , (4i) 

where Iljj is the pressure-strain correlation tensor. 
Following ?, in the case of isotropic turbulence, and with 
their notations, we have 



3 ,2 

n r = 2^ 2U r ^o. 



(42) 



We can derive a formal expression for the turbulent 
diffusion of momentum, using the classical definition , 



u„u A 



Vt = 



CLu 2 



rctO rctO ' 



(43) 
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For steady turbulence, in the limit where molecular 
viscosity in negligible compared to turbulent transport, 
from equation (|41|) , 

2 




3 a. ft 



(44) 



? derived an expression for the turbulent viscosity from 
? and ? experimental data, 



(3Re* 



(45) 



(46) 



v t |r 3 <9 r f2| 

P 

v v 

where (3 was found to be independent from the geome- 
try and the background flow (see also ? and ? for other 
arguments in favor of this prescription). Identifying this 
expression with Eq. 14411 . one finds, 

, / C^-(3/2)C 2 V ( A" " 

? also identified the Re* parameter (named "gradient 
Reynolds number") as the control parameter for stabil- 
ity. This derivation was acheived by identifying the values 
of the Reynolds number for which the measured torque 
from laboratory experiments starts to differ from its the- 
oretical value for the laminar Couette flow. For small gap 
Couette- Taylor apparatus, the classical Reynolds number 
defined as Re = AVLARR/v (where AR is the gap between 
the cylinders, AQ their differential angular velocity, and 
R the mean radius) was identified as the relevant control 
parameter. For gaps bigger than AR/R > 1/20 the criti- 
cal Reynolds number behave like (AR/R) 2 , which means 
that Re* takes a constant value. To be more proper, this 
is a measure of the efficiency of the turbulent transport. 
The bifurcation identified from the torque data is then 
not directly linked to the stability of the flow, but is the 
threshold for which turbulent transport becomes appre- 
ciably more efficient than molecular viscosity. It should 
be read as the following condition 

— » 1 Re* » 4 (47) 
v p 

This is consistent with the numerical values from ? 
(Re* = 1.5 • ICT 5 and = 1.5 • 1(T 6 ) 

5. Geometrical constraints 

5.1. General case 

In the case of a bounded flow, keeping up the assumption 
of isotropic turbulence, we have to assume that A will be 
of order of a macro-scale d. The stability criteria for a vis- 
cous flow is modified as A is constrained to smaller values 
than for a free shear flow (therefore being closer to the 
inviscid values Ro + and Ro~ , according to Eqs l|3"?|) and 
1)38(1 .) while for the inviscid case it remains unchanged. 
The radial transport will also be modified, as we have, 
according to eq. ifSjl 

2 r<? r ft| 



5.2. Narrow gap Couette-Taylor experiment 

This expression can be approximated in the case of 
narrow-gap Couette-Taylor experiment by taking for d the 
distance between the two cylinders AR. 



v t Ail ■ AR ■ R 

— oc = Re, 

v v 



(49) 



cx df 



(48) 



The bifurcation threshold identified on torque mea- 
surement aforementioned will then scale as the classical 
Reynolds number, in agreement with the results from ?. 



6. Discussion 

We have shown that a simple model of turbulent fluctu- 
ations energetics can explain laboratory experiment re- 
sults on differentially rotating flow. The equations of mo- 
tion show that shear instabilities can develop in rotating 
flows, as long as the non-linearities are taken into account. 
Making reasonable assumptions about the necessary prop- 
erties of the velocity correlation products, we derived for- 
mal expressions for the control parameters for stability. 
Comparison with laboratory results confirms that the rel- 
evant parameter is the Rossby number. In the case of a 
viscous flow, the critical Rossby numbers exhibit a ft -1 
correction, therefore converging toward the same constant 
critical Rossby numbers as in the inviscid case, for large ro- 
tation rates. The turbulent viscosity inferred is consistent 
with the prescription previously proposed by ?. The criti- 
cal amplitude of perturbations are also consistent with lin- 
ear stability theory results in the inviscid limit. We would 
like to stress the importance of the pressure fluctuations in 
the growth and self-sustaining of this kind of turbulence. 
Even though the time average of the related terms vanish 
from the velocity fluctuation equations, they can be seen 
as the main engine for isotropizing the turbulence, hence 
disappearing in time average when isotropy is acheived. 

Complementary work is needed, in the laboratory or 
numerically, in order to seek whether the two particular 
values or the critical Rossby number matching the ex- 
perimental results on the Couette-Taylor experiment are 
genuine to that system or are of more "universal" value. 
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Appendix A: Turbulent fluctuations equations 

The turbulent transport equations are given here in cylin- 
drical coordinates (r, <f>, z). Viscous diffusion is omitted 
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o t u r + u r d r u r H a^,u r + u z d z u r - 2 — u r u, = 



— 2 [ u r o r u r + u r u^ I ] + u r u z o z u r 



- 2u' r y-d r (ru' r 2 ) + -84,(11,1/4,) - — + d z (u' r u' z ) 



+ 2-p'd r u' r — 2-d r (p'u' r 



- v \2(d r u r ') 2 + (<9 r u/) 2 + (d r u z ') 2 



/ \ 2 



+ I -dtfUr — \ + {d z U r 'Y 



(A.l) 



dtu'^ + UrdrU^ + —d^u'^ +u z d z u' (t) -Y2 — U r U^ = 

- 2 (u'rU^drUj + U 2 + 7") + U '<t, U 'z d ^ 

(1,2 , u r M rf> / \ 



<9 t ?4 2 + u r d r u' 2 + -^-d^u' 2 + u z d z u' 2 — 

,, 1 - — r a — . — — rdcj,u^ — T2„ — 
2 I u r u z o r u z + u z u, h m z o z u z 



2u z (^-d r {ru' r u' z ) + + d z u' 2 ^j 

2-p'd z u' z - 2-d z {p'u' z ) 
P P 

u(2(d z u' z f + -^(d^f + (d r u' z f 



+{d z u r f + (d z u'^j 



(A.3) 



